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We give a method of describing thermodynamical transport phe- 
nomena, based on a quantum scattering theoretical approach. We 
0^ \ consider a quantum system of particles connected to thermodynam- 

ical reservoirs by leads. The effects of the reservoirs are imposed as 
J> ■ an asymptotic condition at the end of the leads. We derive an ex- 

, pression for a current of a conserved quantity, which is independent 

\ of the details of the Hamiltonian operator. The Landauer formula 

' and its generalizations are derived from this method. 

"j^ i (1. Introduction } Statistical mechanical description of thermo- 

^ ' dynamical responses has been one of the important subjects of nonequilibrium 

' O ' statistical mechanics. Some methods have been proposed for this purpose (for 

O ■ example, see Ref. [1]). Of special interest is a method pioneered by Landauer 



[2], who heuristically derived an expression for an electric current, employing a 
scattering theoretical approach. His method was generalized to some cases; for 
example, a multichannel case [3-5], a case of a finite temperature [4,6], a case 
of a heat current [6,7] and a case of an inelastic scattering process caused by 
a random potential [8]. These methods describe linear responses of a system 
to thermodynamical gaps of reservoirs which induce Fermi distributions in the 
system. However, these have not explicitly treated an effect of inelastic scatter- 
ing processes caused by scatterers, and applications of these methods have been 
mainly restricted to mesoscopic phenomena. Moreover, there have not existed 
a unified statistical mechanical derivation of all the generalizations of Landauer 
formula [9]. 

The purpose of the present Letter is to give a statistical mechanical method 
for descriptions of thermodynamical responses, based on a quantum scattering 
theoretical approach. We show that Landauer formula and its generalizations 
are derived by this method. This method covers all the cases which have been 
contained in the generalizations of Landauer formula. Moreover, it can be ap- 



plied to some new cases; nonlinear responses to thermodynamical gaps of more 
general reservoirs inducing non-Fermi distributions in a system, inelastic pro- 
cesses caused by scatterers, and currents other than an electric or a heat current, 
etc. So, this method can give new generalizations of Landauer formula. 

( 2. Set-up } We consider a quantum system of particles in a three- 
dimensional region fl. The system consists of two kinds of particles, which 
we call 'transport particles' and 'scatterers'. The transport particles are in a 
scattering state, and the scatterers are in a bound state. The region consists 
of a finite region fio and N semi-infinite columned regions ^j, j = 1, 2, • • • , A^. 
The semi-infinite columned region Q,j connects the region to infinity. We 
call the columned region the 'lead'. 

For simplicity, we treat a system of only two particles; one transport particle 
and one scatterer. We assume the Hamiltonian operator H of this system to be 
of the form 

where c is the velocity of light, m and M are the masses of the transport particle 
and the scatterer, respectively, q and Q are the charges of the transport particle 
and the scatterer, respectively, x and X are the coordinate operators of the 
transport particle and the scatterer, respectively, p and P are the momentum 
operators of the transport particle and the scatterer, respectively, A{x) and 
A{X) are the vector potential operators acting on the transport particle and 
the scatterer, respectively, U {x, X) is the potential operator of the transport 
particle and the scatterer. Here the square of a vector means the inner product 
of the vector with itself. 

The state of this system at time t is described by a density operator p(t) 
which obeys the Liouville-von Neuman equation 

^n^ = [H,m (2) 

where 27r?i is the Planck constant. 

We introduce | x,X ) as the eigenstate of the operators x and X with 
eigenvalues x and X, respectively. We introduce the unit vectors e^^\ k = 1, 2, 3 
as a basis of such that e^^ is parallel to the j-th columned region and is 



pointing to the finite region Hq. We define x'j^' by a;^^' = e'j^' ■ x, and introduce 
a;)f ' as a eigenvalue of ai^f' {j = 1, 2, ■ • • , TV, /c = 1,2,3). We assume the functions 
A{x) and i7(a;, X) to have the asymptotic forms satisfying 



A{x) '~ A"'°°'(x^^\4^'), inxeflj, (3) 



U) _ 

U{x,X)''' ^ ^ U^^'"'\xf,xf,X), mxe9.j (4) 

where ) is a function oi and and U^^'°°\x^i\xf ,X) is 

a function of ' , * and X only. We consider the operator H'-'-'^^ defined by 



The eigenstate | 4'Jf^~' ) of the operator ~) can be represented as 



I ^n' ) = \^V)®\ ^tn' )■ (6) 

Here, | 0^^' ) is the eigenstate of the operator e^^' • p with the eigenvalue %k and 
have an orthonormality 



{<t>f \<l>'l!)=2n5{k-k'). (7) 

And I </'$j'„°°' ) is introduced as the eigenstate of the operator which is 

defined by ( (t)f \ | ) = nf^\ (t)f \ (f)'^^ ). We introduce E^^^] as the 

eigenvalue of the operator H'^^-°^'> corresponding to the eigenstate | $5j^„°°' }. And 
I X2\x''^\X ) is introduced as the eigenstate of the operators X2', 2:3' and X 
with eigenvalues Xj', 2:3' and X, respectively. We assume an orthonormality 

/ da;«' dx'f [ dX{ ^Jf;r' | x'f , xf ,X){ xf , xf ,X\^'i^r') = S^n' (8) 
JSj Jn 

of the states {| i^Jf^^' )}«, where Sj represents the projection of the cross section 
of the j-th columned region onto the 3 plane. 



We assume that there exists an eigenstate | ^j.^^ ) of the Hamiltonian op- 
erator H with the eigenvalue E^^^^. Here, the eigenstate | "^^j^l^ ) satisfies the 
following asymptotic condition: 



k n 
kk' <0 



in cc G f2j 



(9) 



k' n' 
kk'<0 



A') 



in a; G , I ^ j 



where r'fl , , , and tj'"/' , , , are constants determined by the fact that I ^'Sf ' ) 
is an eigenstate of the operator H . It should be noted that the eigenstate | ^'Jf^'^ ) 
may not be determined uniquely by the asymptotic condition (|^) . So we should 
introduce a new suffix in order to distinguish the eigenstates having the same 
asymptotic form. But in order to avoid a complicated notation we do not write 
such a suffix explicitly, but distinguish such degenerate states by different values 
of n from now on. 

We assume that the condition 

lim ( a;, X I p(t) | A" ) = (10) 

\X I — >oc 

is satisfied for an arbitrary state \ X )■ The condition (|l^ may be satisfied by 
the assumption that the scatterer is in a bound state. The transport particle is 
assumed to be in a scattering state which satisfies the asymptotic condition 



N 



J2 E F,ik,n){x,X\^^:^){^^:^\x',X') (H) 



kn 
{k>0) 



where Fj{k, n) is a positive function of k and n. We interpret that the function 
Fj{k,n) represents a distribution induced in the system by an interaction with 
the reservoir at the end of the j-th lead. 

( 3. Current of a Conserved Quantity ) We consider a conserved 
quantity. This quantity corresponds to a Hermitian operator G, which satisfies 
the relation [H, G] — 0. 

We introduce the symmetrized product X -ky oi two arbitrary operators X 
and y as 



Xi.y = {xy + {xy)''}/2 (12) 

where (Xy)^ means the Hermitian conjugate operator of the operator Xy. By 
Eqs. (H), (|l^), ( p^ ) and Gauss' divergence theorem we obtain an equation of 
continuity 



^JjX{x,X I {p{t)*G)\x,X) 

+ [ dX{x,X I {{p{t)*G)*v) \x,X)=0 (13) 
ox 

where i) is defined hy v = — [H,x]/{ih). By Eq. ( p^ we interpret the quantity 
J^dX{ x,X I {{p{t) -k G) -k i)) \ x, X ) as a. current density of the conserved 
quantity due to the transport particle at time t. 

We consider the current Jj in the j-th lead, defined by 



Jj — lim 



/ dx<^'d4^"' e<^'- / dX{x,X \{{p{t)i.G)i<v)\x,X).{U) 
) Jsj Jn 



We assume that the eigenstate | ^'j^;'^'^ ) is also the eigenstate of the operator G 
[14]. And we introduce Gj^?^ as the eigenvalue of the operator G corresponding 
to the eigenstate | "^''-^^ ). From Eqs. (Ill]) and (III) we derive 



kn k' n' 

(fc>0) / fc'>0 



U) |2 O'^-k-n- 



— i-'r,' kn 



l — l kn k'n' 

{l^j) (fc>0) f k'>0 



where we used the foUowing fact [15]: 



If El'^ = and kk' > 0, then 

dx'^''dxf I dX 
Jo. 

\{{x,X\{e^^ ■v)\^tn'){^t::^\^^X) 

T^^^n'nSk'k, inxerij. (16) 
h ok 

By the equation of continuity (|l^), the current of the quantity G which flows 
through any cross section of the j-th columned region takes the value Jj in any 
steady state. It is important to note that the effect of interference among the 
incident wave, its reflective waves and its transmitted waves does not appear in 
the quantity Jj , because of Eq. (^6|) . Eq. ( p^ is one of the main results in the 
present Letter. 

( 4. Derivation of Landauer Formula and its Generalizations ) 

As a special case, we consider a case satisfying the following conditions: (A) The 
scattered wave has the same wave number as its incident wave, that is, when 
^[1'.', ^ 0, kk' > and E'^,,^, = in n', fc, n on T,, or tj':;*, ^ 0, 

kk' > and E'^^f^,^, = Ef^^ in k, n on and in n' on !Fi, the relation k' = k is 
satisfied, where J-j is the support of the function Fj{k,n). It implies that the 
scattering processes are almost elastic. (B) The distribution function Fj{k, n) is 
a function only of the energy eigenvalue E^^ . So we put a function Fj {E^^ ) of 
^kn instead of Fj{k, n). (C) The operator G is a function only of the operator 
H. So we put an operator G{H), which is a function of H, instead of G. (D) The 
coefficients r'^! ^, and t''^^^\i.^ are dependent on the suffixes k and n on J-j 
only through the energy eigenstate E^^. So we put the coefficients r'^i ^/-^{E^^) 
and i[fc/'„')(£'{,^^) instead of the coefficients ^, and iJfc/'„/.fc respectively. 
(E) The function E^^.^ of the variable k projects the domain (0, oo) of the variable 
k to the domain Aj of the variable e. The domain Aj is independent of the value 



of the suffix n. (F) The suffix n takes n of values on Tj. Under the conditions 
(A)-(F), Eq. ® becomes 




de F,(£)G(e)(l-i?,(£)) 

A, 



J2 ^ d£^i(£)G(e)r,,(e)|. (17) 



1=1 



Sfc'n' I ^[-'k' „')(^) respectively where the sums over the suffixes k' and n' 
are taken only over values satisfying the condition E''^\, , = e and E^^,, , — e, 



Here i?,(e) and T,,{e) are defined by Rj{e) = J^k'n' I fc',„')(e) I' and T,,(£) 

-k'n' —k'n' 

respectively. And we used the fact that the sum over the suffix fc in fc > 
corresponds to the integral over the suffix fc in fc > multiplied the factor 
l/(27r). 

From Eq. ( p^ ) we can derive Landauer formula and its generalizations which 
have already been proposed. For example, iiG{H) = q, Fj{e) = limT-to{exp{(e— 
IJj) / (ksT)} + 1}^^ = 9{e — jij) (where fc^, T and are positive constants), 
I (/ij — fl)/p, \ << 1 (where fi = {fii + fJ,2 + ■ ■ ■ + ^J'N)/N), n = 2 and Aj — (0, oo). 



then Eq. (17) become 



N 

nil ^ — ' 

1=1 

in the first approximation of {^j — p.) /ft- Eq. ( p^ ) is one of the generalizations 
of Landauer formula, and was proposed by M. Biittiker [5]. 

( 5. Conclusion and Remarks ) In the present Letter, we have 
described thermodynamical nonlinear responses of a quantum system to ther- 
modynamical gaps of reservoirs. Here, the thermodynamical reservoirs were 
connected to the system with leads, and effects of the reservoirs were intro- 
duced as an asymptotic condition at the ends of the leads. Thermodynamical 
gaps of the reservoirs cause thermodynamical transport phenomena. We derived 
an expression of a current of a conserved quantity, which is independent of the 
details of the Hamiltonian operator. For example, this method can describe 



an energy current which is caused inside a system connected to heat reservoirs 
having different temperatures. 

This method leads to Landauer formula and its generalizations which have 
already been proposed. Moreover, we can also give further generalizations of 
Landauer formula by this method. For example, in this method we can deal 
with nonlinear responses, reservoirs inducing non-Fermi distribution in a system, 
inelastic processes caused by scatterers and currents other than an electric or 
a heat current, etc., which have not been treated earlier. The method in the 
present Letter can be generalized to the case of many transport particles and 
many scatterers [16]. (However the two- particles system, which we discussed in 
the present Letter, can be interpreted as a mean field approximation for a system 
consisting of many transport particles and many scatterers.) This method can 
also be generalized to the case where interactions of the particles are given by 
a complex potential only in a finite region. 

We interpreted that the function Fj{k, n) of k and n represents a distribution 
induced in the system by an interaction with thermodynamical reservoir at the 
end of the j-th lead, because it represents a distribution of the incident wave at 
the end of the j-th lead, and we can use an equilibrium distribution function as 
the stationary distribution function Fj{k,n). But these reservoirs do not cause 
a dissipation of the system in this method. And as the distribution function 
Fj{k, n) we can select a distribution function which is not an equilibrium distri- 
bution function. In this sense, we can describe non-thermodynamical transport 
phenomena by this method. 

One may notice that this method does not treat distributions of thermody- 
namic quantities inside a system, because we do not assume the local equilib- 
rium assumption, etc. in this method. For example, this method does not treat 
temperature distribution inside a system connected to heat reservoirs having 
different temperatures. To treat such a distribution by a generalization of this 
method remains as a future problem. 

( Acknowledgements ) I wish to express my gratitude to K. Kitahara 
and T. Hara for a careful reading of the manuscript and for valuable comments. 
I acknowledge S. Komiyama for a stimulating lecture on Landauer formula. 



[I] D. N. Zubarcv, Nonequiliblium Statistical Mechanics, Consultant Bureau, 
New York (1971) 

[2] R. Landauer, Philos. Mag. 21 (1970) 863 
[3] M. Y. Azbel, J. Phys. C14 (1981) L225 

[4] M. Buttiker, Y. Imry, R. Landauer and S. Pinhas, Phys. Rev. B31 (1985) 
6207 

[5] M. Biittiker, Phys. Rev. Lett. 57 (1986) 1761 

[6] H. L. Engquist and P. W. Anderson, Phys. Rev. B24 (1981) 1151 

[7] U. Sivan and Y. Imry, Phys. Rev. B33 (1986) 551 

[8] Y. Gefen and G. Shon, Phys. Rev. B30 (1984) 7323 

[9] There is a statistical mechanical derivation of a generalization of Landauer 
formula for an electric current in a restricted situation where the system consists 
of non-interacting electrons, reservoirs connected to the system do not have 
temperature gaps, and the system causes only elastic scattering processes [10- 
13]. But it should be noted that the methods used in Refs. [10-13] are not same 
with our method. Refs. [10-12] used the linear response theory to an external 
disturbance derived from the potential energy, and Ref. [13] used a different 
asymptotic condition with our method. 

[10] E. N. Economou and C. M. Soukoulis, Phys. Rev. Lett. 46 (1981) 618 

[II] D. S. Fisher and P. A. Lee, Phys. Rev. B23 (1981) 6851 

[12] H. U. Baranger and A. D. Stone, Phys. Rev. B40 (1989) 8169 

[13] S. Komiyama and H. Hirai, Phys. Rev. B54 (1996) 2067 

[14] The exchangeable operators H and G have a common eigenstate, and the 

state I '^''i^D is a eigenstate of the operator H. But the common eigenstate of 

operators H and G may not have the asymptotic form (|^). So, this is a new 

assumption. 

[15] The outline of the derivation of Eq. ( p6| ) is as follows. We first note that 
the divergence of the integrated function of ' and 0:3 ' in the left hand side of 
Eq. (|l6|) must be zero when £[^^1 = E'^,\^,. So the left hand side of Eq. (^ must 
be zero except for the case fc = fc' in | a;^;'' [^ 00, x G ilj when E'j^^ = ^k''n'- 
Besides, we have {x,X\ (e<^' • *) | <^f^'> ) {x,X \ ^{d'H'-^'°°^ /dk) \ $<5„°°' ) 
as I cx) in a? G ilj. By these and Eq. (||), we obtain Eq. (p^). 

[16] For such a generalization we must deal with energy eigenstates of the system 
consisting of the many transport particles and many scatterers, which must be 
required to satisfy the Pauli principle. 



